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ABSTRACT 



This study was conducted at the Naval Postgraduate School 
to investigate the formation of particle flow cells at mini- 
mum fluidization gas flow rates in a rectangular gas-fluidized 
bed. The primary objective is to determine if multiple cell 
formation occurs as the bed width to height ratio increases. 
The settled bed height is maintained at a constant level for 
all runs. The secondary objective is to determine if the 
presence of external plate heaters at the end walls of the 
particle bed will alter the particle flow cell formation 
patterns when energized and to what extent they affect the 
fluid mixing (heat transfer) flow patterns within the par- 
ticle flow. Pressure measurements across the bed are used 
to determine minimum fluidization points and thermocouple 
probe measurements determined the thermal gradients within 
the bed for each bed configuration. Visual observations on 
the external bed movement are included in the determination 
of the results of this experimental study. 
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NOMENCLATURE 



English Letter Symbols 

a Constant, dependent on tube orientation - Reference 30. 

A Cross-sectional area of bed container, 

c 

A^ Contact area of packet and surface - Figure 18. 

B Shape factor - Reference 12. 

3i Biot number. 

C m Specific heat of packet or quiescent bed. 

C ^ Specific heat of fluid. 

Cp S Specific heat of solid particle. 

C R Constant, correction for location of immersed body - 

Reference 34. 

C s Specific heat of solid particle. 

dp Diameter of sphere of same specific area as that 

of the bed particle. 

D Diameter of the fluidized bed. 

D' Effective diameter of free area across the bed. 

D H Diameter of an immersed heater. 

Particle diameter of sieved fraction. 

0^ Particle diameter. 

D Average particle diameter. 

IT 

D Distance between two successive layers - Figure 16. 

D^ Diameter of immersed tube. 

D£q Diameter of immersed tube being equal to 20 mm. 

g Acceleration of gravity. 
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G 

G' 

Ga. 



G 

e 




G v 

h 



h o 

h w 

^wiL 

h 

wL 

n wmax 




Fluid mass velocity. 

Fluid mass velocity based on voidage area. 

Galileo number (Archimedes number) . 

Hypothetical fluid mass velocity required to 
merely expand a bed of particles. 

Fluid mass velocity through a fluidizing bed. 

Fluid mass velocity for minimum fluidization. 

Optimum fluid mass velocity giving a maximum 
value of the heat transfer coefficient. 

Fluid mass velocity based on kinematic viscosity. 

Heat transfer coefficient for the contact 
resistance at the wall. 

Empirical constant - Reference 44. 

Heat transfer coefficient between the fluidized 
bed and the surface. 

Instantaneous local heat transfer coefficient. 

Local average heat transfer coefficient. 

Maximum heat transfer coefficient between 
fluidized bed and surface. 

Superficial bubble-to-cloud heat transfer coefficient 
Fluid thermal conductivity. 

Thermal conductivity of packet or quiescent bed. 
Thermal conductivity of particle. 

Thermal conductivity of sphere. 

Distance from bed top in packet theory - Reference 47 
Fluidized bed height. 

Expanded bed height. 

Length of heat exchange surface. 
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L ^ Bed height at point of minimum fluidization, 
m Slope - Figure 12. 

m" Slope - Figure 15. 

N Fluidization number. 

Nu Nusselt number. 

Nu Maximum value of Nusselt number. 

IUclX 

p Empirical constant - Reference 44. 

P Center-to-center distance (pitch) of adjacent tubes. 

Pr Prandtl number, 

q Heat flux. 

q r Heat transferred into the core portion of the bed by 

an interchange of particles - Figure 23. 

q w Heat flow rate from the wall into a packet. 

q . Instantaneous heat transfer rate. 

W1 

q Heat taken by particles flowing parallel to the 

surface in a second zone of emulsion of thickness 
- Figure 23. 

r Fluctuation ratio - highest to lowest bed level in a 

gas-fluidized charge. 

r S p Radius of sphere. 

R Bed expansion ratio. 

Re Reynolds number. 

Re^ Reynolds number based on particle free fall velocity. 

Re m Modified Reynolds number - Reference 27. 

Re m £ Reynolds number at minimum fluidization. 

Re Optimum value of Reynolds number giving a maximum 

^ heat transfer coefficient. 

Re^ Reynolds number based on bed diameter. 
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u d 
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u mf 

u o 
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p 

U Q 
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s 

u t 

V 
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W 

x 
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1 

X 

y 



Average replacement of packets at wall by side mixing. 
Stanton number. 

Time . 

Bulk bed constant temperature. 

Fluid film temperature. 

Gas temperature into bed. 

Gas temperature out of bed. 

Temperature of wall or heater. 

Initial fluid velocity with no particle motion. 

Fluid velocity required for initial expansion. 

Fluid velocity required for minimum fluidization. 
Superficial velocity at fluidization conditions. 
Particle velocity. 

Fluid velocity during fluidization. 

Particle velocity of Packet Model. 

Terminal fluid velocity. 

Volume displaced by particles. 

Bed width. 

Axis direction - Figures 21 and 22. 

Weight fraction of sieved particles for . 

Packet width - Figures 21 and 22. 

Horizontal pitch of tube bundle - Reference 47. 
Vertical pitch of tube bundle - Reference 47. 



Greek Letter Symbols 
a Constant - Reference 9. 
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Constant - Reference 9 . 



Specific weight of fluid. 

Boundary layer thickness at point between two particle 
layers - Figure 16. 

Average boundary layer thickness at point between two 
points - Figure 16. 

Thickness of emulsion layer. 

Thickness of gas layer. 

Minimum distance between surfaces of adjacent tubes 
(gap) . 

Pressure drop across bed at u^. 

Pressure drop across bed at u £ . 

Pressure drop across bed at u m f. 

Ideal pressure drop across the bed calculated from the 
buoyant weight of the bed. 

Void fraction at fluidization conditions. 

Initial particle voidage as filled. 

Particle voidage during expansion. 

Particle voidage at minimum fluidization. 

Void fraction at fluidization conditions. 

Fluidization efficiency. 

Non-dimensional absolute bed temperature. 

Fluid viscosity. 

Fluid kinematic viscosity at minimum fluidization. 
Constant . 

Fluid density. 

Density of packet. 

Fluid density at minimum fluidization. 



p Particle density, 

s 

P s k Nonfluidized bed density. 

<() Sphericity of solid particle. 
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INTRODUCTION 



I . 

If one considers a fluid flowing through a bed of solid 
particles, several forms of today's fluidized bed technology 
are created. When the particles are stationary, the process 
is called flow through a porous medium or fixed bed fluidiza- 
tion. If the fluid velocity is such that the particles 
become entrained by the fluid and are removed from the bed, 
the process is then called pneumatic conveying or entrained 
flow fluidization. At some point between these two conditions, 
the particles in mixture with the fluid are transformed into 
a fluid-like state causing continuous interaction between the 
two mediums. The bed particles in this state do not, however, 
leave the bed container. This state is called fluidized bed 
fluidization . 

The technique of fluidization (all forms) is used for 
solid-fluid interaction processes v/henever high rates of heat 
and mass transfer between the two constituents is required or 
high rates of particle transport to and from the bed is 
necessary. The advantages in the use of fluidization can be 
listed as follows: [6,12] 1 

(a) The formation of an extremely large area of contact 
between the particles and fluid. This very large surface 

lumbers in brackets indicate references. 
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area permits the achievement of high overall rates of heat 
and mass transfer between the particles and the fluid. 

(b) The comparative ease with which the fluidized 
particles can be handled. 

(c) The reduction of temperature gradients to negligible 
proportions throughout the bulk of the bed as a result of the 
high degree of particle mixing that can occur in gas-fluidized 
systems, i.e. a very high effective internal thermal conduc- 
tivity . 

(d) The occurrence of high rates of heat transfer between 
the fluidized particles and an immersed surface. 

(e) The value of the high thermal inertia of the solids 
in heat transfer operations within the bed as a whole. 

(f) The face that the gas-fluidized systems approximate 
a fluid which has a low vapor pressure even at very high tem- 
peratures - as a thermodynamic fluid it lies roughly between 
gases and liquids. 

Inherent with the fluidization techniques are also some 
limitations: [6,12] 

(a) The process demands the expenditure of power for 
fluidization. 

(b) Some solids cannot be fluidized because of either 
an excessive tendency for particle attrition or agglomera- 
tion to occur. The first will produce too many fines by 
degradation of the particles which will be lost from the bed. 
The latter will give "agglomerated” particles too large to 
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fluidize. Also, there are those materials of such irregular 
shape that they can only be handled by spouting bed techniques. 

(c) Operating rates are limited to within the range over 
which the bed may be fluidized. If the fluid flow is too low, 
the bed settles; if, on the other hand, the rate is too high, 
there is an excessive loss of particles from the bed by carry- 
over . 

(d) There is usually a limit to the size of the particles 
which can be fluidized. This limit is approximately 6.5 mm 
(0.25 in) down to a few microns (0.001 in). It is particularly 
dependent on the fluid viscosity and particle/fluid density 
ratio as well as the fluid flow rate. 

(e) True counter flow operation is unobtainable because 
of the high degree of mixing that usually occurs within a 
fluidized bed. 

(f) The general absence of thermal gradients throughout 
a well mixed fluidized bed is a disadvantage in those opera- 
tions where thermal gradients are required. 

(g) The dynamics of fluidized systems are not sufficiently 
well understood to permit confident scale-up to large scale 
units from small scale test data when reaction conditions are 
critical . 

Up until this time the term fluid has been used in describ- 
ing the transport constituent in fluidization. To be more 
specific there are two fluids which have received widespread 
use in fluidized bed technology, air (gas) and water (liquid) . 
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The basic differences are in the dynamics and the mass and 
heat transfer of the solid particles suspended in the gas- 
fluidized or liquid-fluidized bed. This is primarily a result 
of the difference in the particle to fluid density ratios. In 
the gas-fluidized bed the ratio will usually be on the order of 
a thousand while in a liquid-fluidized bed the particles are 
only several times heavier than the fluid. 

Additionally, the fluid flow characteristics will affect 
the overall mass and heat transfer bed characteristics. In 
the liquid-fluidized system, the fluidization velocities are 
usually of such a magnitude that the fluid inertial terms are 
considered insignificant. This is not the case with a gas- 
fluidized system, especially when large particle sizes are 
used. As the velocities increase the inertial fluid losses 
and form drag become important in determining the bed charac- 
teristics . 

With the smooth manner in which liquid-fluidized beds 
operate, i.e. no gross instabilities in flow, the bed is 
sometimes referred to as a homogeneously fluidized bed or 
a particulately fluidized bed. The gas-fluidized beds on 
the other hand usually have unstable flow above some minimum 
fluidization condition, characterized by violent agitation 
of the particles accompanied by bubbling and channeling of 
the gas. As a result, the gas-fluidized bed is sometimes 
called an aggregate or bubbling fluidized bed. 
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The bulk of recent research has been in the area of 



gas-fluidized beds vice liquid- fluidized beds . Whether this 
trend will shift remains to be seen. Nevertheless, this 
investigation will be in the gas-fluidized field and all 
further discussions will be based on available gas-solid 
research data. 

In a gas-fluidized bed the particles are supported above 
a gas distributor plate which permits the gas to enter and 
flow upward through the bed. As the gas flow increases, 
there will be a pressure drop across the bed which will 
eventually reach a point when it just balances the downward 
gravitational force of the particle bed. When this point is 
reached, the bed will expand and the particles separate. 

This is referred to as the minimum fluidization point. Fur- 
ther gas flow rate increases will result in bubble formation 
which will rise through the bed and burst at the surface not 
unlike a vapor bubble in a boiling liquid. 

The separation of the particles creates a large particle 
surface in direct contact with the gas and the bubble flow 
induces vigorous and thorough mixing of the particles. This 
high gas-particle contact area results in the high heat 
transfer rates which have been experienced. Because the 
particle flow takes place in the fluidizing gas, a secondary 
fluid-mixing phenomenon will be induced through the particle 
flow pattern. The extent of this secondary mixing and its 
effect on heat transfer will depend on whether the particles 
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are dispersed in the bed more or less as individual particles 
or clumps of particles. As an individual particle, disper- 
sion by the fluid mixing would be relatively small. As an 
agglomerate of particles, however, the result is a dragging 
of gas pockets through the bed resulting in perturbations to 
the particle flow patterns. Thus, the extent of fluid mixing 
is related to the mode and quality of fluidization. This is 
an area that has been scarcely investigated. Two properties, 
though, that have been established as a result of fluid 
mixing investigations are that: 

(a) Owing to fluid back-mixing, temperature driving 
forces are reduced. 

(b) Analogously to heat transfer, driving forces are 
also reduced with mass transfer. 

Based on previous research the normal particle flow 
pattern in small diameter cylindrical gas-fluidized beds is 
as shown in Figure 1. For gas flow rates which are signif- 
icantly greater than the minimum fluidization flow rate the 
solids tend to move downward at the bed walls and the passage 
of the individual particle is readily observed. The particles 
will move down the wall and then submerge into the body of 
the bed. Around this entry point other particles will emerge 
from the bed and travel down the wall until they also merge 
back into the bed. This process continues until the bottom 
of the bed is approached. At some point near the bottom of 
the bed an inactive or dead region will form in the particles. 
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The size of this region will depend on the height of the bed, 
fluid flow rate and bed container geometry. In deep beds 
the flow regions are well defined while in very shallow beds 
the downward flow pattern is much less clearly defined with 
the bed appearing to be more turbulent and in a homogenous 
state of agitation. 

In the literature reviewed, when the term shallow bed was 
used it generally was referring to a particle bed whose depth 
was on the order of one and a half feet or less depending on 
the fluid inlet cross-sectional area. In general, the major 
features distinguishing shallow beds from deep beds can be 
described as: [55] 

(a) Violent coalescence of bubbles taking place through 
the majority of the shallow bed depth. 

(b) The jetting region next to the distributor plate 
occupying a significant proportion of the shallow bed depth. 

(c) The thermal "entry region" occupying a significant 
proportion of the shallow bed depth. 

These features would suggest that the particle motion 
within the shallow bed may well be significantly greater 
than in a deep bed and that the bubble growth will largely 
be completed within a very small proportion of the bed depth. 
Thus the residence times of the particles at the heat trans- 
fer surface would be considerably shorter than in a deep bed. 

In the flow description of deep beds, the bed particles 
will move upward in the center of the bed either in a spouting 
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type pattern or in a multi-channel pattern. It should be 
noted that if the bed vessel i.s operated in a configuration 
such that its containment walls are not parallel to the 
fluid distributor plate inlet gas flow direction, then the 
particle flow pattern will be adversely affected. Only a 
few degrees of deviation from this parallel configuration 
will lead to particle flow pattern short-circuiting. If the 
operating gas velocities are further increased while this 
configuration is present, then slugging will initialize. 

As was previously stated, fluid mixing in gas-fluidized 
beds will occur as a result of the particle mixing at the 
onset of fluidization. Unfortunately, the gas does not 
circulate in the smooth manner shown in Figure 2. Instead, 
there will be many eddies within the central core of the 
particle flow pattern as well as horizontal components 
radiating out from the core. Additionally, part of the gas 
may rise rapidly in the form of bubbles, rising ahead of the 
slower moving interstitial gas components. Therefore, 

Figure 2 should only be viewed as an ideal flow pattern at 
minimum fluidization for the bed. 

To date, particle flow patterns with the resulting par- 
ticle and fluid mixing criteria have been explored only in 
equipment of relatively small diameter (less than one foot) . 
The experimental apparatuses that were used were circular 
in nature with either no external heat source or a submerged 
cylindrical (tube) heat source. The data obtained was 
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Figure 1. Observed particle circulation pattern. 
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Figure 2. Greatly simplified gas f 1 out pattern. 



26 



usually a byproduct of the investigation to determine heat 
transfer coefficients for fluidized beds with various heat 
sources and bed geometries . 

This thesis investigation was conducted to determine 
the particle flow patterns developed in a gas-fluidized bed 
of variable rectangular geometry. Secondly, to determine 
what effect an external plate heat source has on the particle 
and fluid flow patterns developed within the rectangular bed. 
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II. OBJECTIVES AND METHOD OF APPROACH 



A. GENERAL STATEMENT OF THE PROBLEM 

The specific problem may be defined as an investigation 
into the formation of particle flow cells at minimum gas 
fluidization flow rates in a rectangular gas-fluidized bed. 
The primary objective is to determine if multiple cell for- 
mation occurs as the bed width to height (W/L e ) ratio 
increases. The secondary objective is to determine if the 
presence of plate heaters at the end walls of the particle 
bed will alter the particle flow cell formation patterns 
when energized and to what extent they affect the fluid 
mixing (heat transfer) flow patterns within the particle 
flow. 

B. METHOD OF APPROACH 

The methods used for this investigation were experimental 
and analytical. Using available research literature a deter- 
mination was made of the minimum fluidizing velocity (u mf ) 
for the particular fluidized bed configuration. This value 
of u ^ was based on recorded and analytical pressure drops 
across the distributor plate and bed in addition to the 
expanded bed heights. Visual confirmation of the minimum 
fluidizing conditions was then compared with the analytical 
estimates of the required pressure drops, bed height and air 
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velocities. Once an initial particle flow cell was estab- 
lished, the bed geometry and/or activation of the plate 
heaters would be used in an attempt to alter its configura- 
tion. A grid of variable depth copper-constantan thermo- 
couples was used to map the thermal distribution within the 
particle flow cell to determine the fluid mixing properties 
of the cells. 

This thesis will not be involved with the heat transfer 
coefficient determination so prevalent in recent experimenta- 
tion; rather it will focus only on how the particle flow cells 
relate to the bed geometry chosen and how the resulting fluid 
mixing affects the heat transfer within the bed. It should 
be noted that direct measurement of the solid particle temper- 
atures is very unreliable and that the temperatures recorded 
will be considered to be that of the air at a particular 
point within the grid of the bed. Because of the difference 
in the particle to fluid density, however, the two will be 
considered to be at equilibrium at the point measured for 
this experiment. 
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III. SUMMARY OF PREVIOUS RESEARCH IN 
GAS-FLUIDIZED BED TECHNOLOGY 



A. GENERAL DESCRIPTION OF GAS-FLUIDIZED BED BEHAVIOR 
1 . The Phenomenon of Fluidization 

When the bed configuration is initially filled with 
particles, an initial voidage (e^) will result between the 
particles. As the fluidizing fluid, a gas in this case, is 
initiated at some superficial velocity (u^) there will be no 
movement of the bed. However, a pressure drop (Ap^) will 
occur and gradually increase at a linear rate. The magnitude 
of this pressure drop will be determined by the gas flow rate 
and the characteristics of the bed. 

Eventually a condition will be reached when the gas 
flow rate (u e ) has increased to a point at which expansion 
of the bed initiates. The voidage will increase to some e Q 
and the pressure drop can be calculated from the equation 

AP e = L e ( 1 " e e } ( p s " p f ) (3.1) 

(Note - The units to be used with all text equations can be 
found in Appendix D.) While the bed is expanding, there 
still will have not been any particle motion taking place. 

As the gas flow rate is increased further, a point 
will next be reached when the voidage becomes and the 

particles are at the point of fluidization onset. The 
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velocity is called the minimum fluidization velocity (u mf ) 
and is critical to most bed behavior calculations. Further 
gas velocity increases (u^) will fluidize the bed causing 
further voidage increases (e^) until a terminal velocity 
(u fc ) is reached as defined by the bed container geometry or 
particle entrainment takes place. By the time u fc has been 
reached the particle motion will be intense although the 
pressure drop across the bed will remain at the level of 
minimum fluidization inception. This condition is referred 
to as dense-pnase turbulent fluidization at its inception 
(Uq) but degenerates to a dilute-phase when the gas rate is 
such that the particles become entrained in the gas. On the 
average the ratio of the terminal velocity (u t ) to the 
minimum fluidization velocity (u m f) is about seventy (70) . 
The terminal velocity is also a function of the particle 
diameter, shape and density as well as the properties of the 
fluidizing fluid. 

Parent et al . [52] has characterized three types of 

behavior with the gas-solid mixtures in the dense-pnase 
fluidization region. Under certain conditions, usually for 
small particle sizes, the gas tends to channel through the 
bed without causing movement of the bed. Under other condi- 
tions, such as large bed heights, large particle sizes and 
high gas flow rates, slugging takes place which is analogous 
to the turbulence observed in a boiling liquid. The third 
type of behavior is described as the mechanically smooth 
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fluidization. In this case the particles move along twisting, 
ever-changing paths. Occasionallly gas bubbles appear in the 
bed but there is no rough bumping. It is this case that the 
pressure drop across the bed is nearly equal to the suspension 
head opposing the flow (drag) . A word of caution is warranted, 
however, in that in research conducted by Lewis et al . [42] 
the conclusion that the pressure drop is essentially equal to 
the weight of the particle bed would appear to be satisfactory 
only for units with small L/D ratios. It appears to be in 
considerable error for large L/D values or for beds having 
large static effects present. 

Here again several researchers differ in their opin- 
ion as to the cause of this phenomenon. While Lewis et al . 

[42] explained it as an indication of frictional drag of the 
solid particles on the walls of the unit, Toomey and Johnstone 
[62] attribute it to kinetic energy losses caused by the colli- 
sion between the particles and Wilhelm and Kwauk [76] con- 
cluded it was caused by forces which are electrostatic and/or 
fluid dynamic in origin. Since this thesis is operating in 
the shallow bed region it is not anticipated that this excess 
pressure drop will be of concern. 

2 . Voidaqe at Minimum Fluidization 

The first bed characteristics that must be determined 
for an understanding of gas-fluidized bed operation is that 
of the bed voidage at minimum fluidization (e^) • By defini- 
tion is that particle separation to which the bed must 
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expand before particle motion, as induced by fluid flow, can 
initialize. There are two methods which can be used to deter- 
mine this factor. 

The first method is based on a collection of data from 
several researchers that has been graphically combined by Leva 
[36] and reproduced here in Figure 3. For a given particle 
diameter (D ) and composition a value for the minimum fluidiza- 

IT 

tion voidage (e m f) tan be obtained. The obvious drawback for 
this method is that it depends solely on the material chosen 
for the particles being one that has had previous experimenta- 
tion accomplished. 

The second method is used with those beds for which 
particle data is not available. It requires that the bed 

be subjected to the fluidizing gas until particle motion 
initializes, at which time the bed height (L^) is recorded. 

The minimum fluidization voidage (e can then be calculated 
from the equation 



V Y* 
p 'f 



e . = 1 - 
mf 



L , A ( p - p- ) 
mf c s f 



(3.2) 



The accuracy of this estimate will depend on the mass charac- 
teristics of the bed particles, i.e. sphericity, whether solid 
or hollow, degree of particle porosity, etc. Particles of 
high porosity are usually unsuitable for fluidized bed use 
as they have a tendency to break up with increased fluid flow. 
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a-So-ft Brick 

b-Rb sorption Carbon 

c-Broken Raschig Rings 

d-Coal and Glass Powder 

e -Car bo rundum 

f-Sand 




Dp (in.) 



Figure 3. Values o-f C vs. particle diameter. 



. 58 ) 



34 



Once a value for the minimum fluidization voidage 



has been estimated, additional bed behavior characteristics 
can be predicted and the bed's suitability for the intended 
process determined. 

3 . Pressure Drop 

Another excellent indication of fluidization behavior 
of the particle material may be obtained from the pressure 
drop-flow relationship of the expanded bed [36,40,64], Ideally, 
the pressure drop across the bed should increase linearly until 
the voidage has reached a maximum value (e ) at which time the 
pressure differential would level off for the remaining phases 
of fluidization. See Figure 4. 

Unfortunately, the idealization never occurs as a 
result of three major factors; namely channeling, slugging 
and distributor design. With channeling a pressure drop 
occurs (Figure 5) when e e is reached. The magnitude of this 
drop will depend on the extent of the channeling phenomenon 
but eventually a minimum drop point is reached and the pres- 
sure differential then increases to some new lower steady 
state value. 

The channeling can be either "through channeling" in 
which it extends completely through the bed or "intermittent 
channeling" where it extends through part of the bed. For- 
tunately, as the gas flow rate is increased the tendency for 
the bed to have channeling decreases. Some factors [38] 
which tend to promote channeling are poor distributor design. 
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Figure 5. Pressure-drop-f 1 oui diagram -For channeling. 
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the shape of the bed particles (less sphericity the greater 
the channeling) , moisture in the bed and the rate of fluid 
flow through beds of small particles. If the bed is moist, 
the particles will tend to clump together, and the gas flow 
will channel between the clumps. At higher flow rates the 
increased agitation of the bed apparently destroys the 
channels as soon as they are formed. 

With slugging a condition exists in which the bubbles 
of gas combine to a size approaching the order of magnitude 
of the diameter of the confining vessel. As a result, the 
particle layers or slugs of granular particles will move 
upward through the bed like a piston until they disintegrate 
at the bed top. When this happens, the pressure drop will 
increase above the predicted steady state value (Ap g ) , see 
Figure 6. If the flow rate is increased further, the pres- 
sure drop rise continues and the bed will slug even more 
severely . 

The third factor, namely distributor design, was 
well investigated by Grohse [27] who stated that for the mul- 
tiorifice plate the point of incipient fluidization is not 
well defined. The gradual pressure-drop increase up to the 
theoretically demanded suggests that the bottom portion of 
the bed is not fluidized to the same extent as the upper 
region. With the screen distributor, a certain amount of 
channeling probably occurred. There is the possibility that 
the holes in the screen were locally blocked; this could 
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Figure S. Pressu re-d rop—f 1 ou diagram for slugging. 
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create a condition favorable for channeling. With the porous 
plate the point of incipient fluidization is well defined 
and there is no tendency for channeling. See Figure 7. 

Since the point of minimum fluidization is the key 
starting point on fluidization studies, it seems extremely 
important to be able to predict that point both theoretically 
and through experimental results. Therefore, the distributor 
design's effect on the pressure drop at u m ^ becomes a very 
important factor. According to Kunii and Levenspiel [34], 
distributors should have a sufficient pressure drop to 
achieve equal flow through all openings. As a result, the 
pressure drop across the distributor plate by direct measure- 
ment should be at least 10% of the pressure drop of the bed 
to ensure effective fluidization. Theoretically the pressure 
drop at the minimum fluidization velocity (Ap^) can be 
found by substituting L mf and for the values of L e and 

e g in equation (3.1) . 

4 . Minimum Fluidization Velocity 

The ability to predict the point of initial fluidi- 
zation reliably is of basic importance in virtually all 
f luidized-process studies and designs. By operating at mass 
velocities which are slightly higher than the minimum fluid- 
ization mass velocity (G^) attrition and elutriation 
(stratification) losses will be minimized. The by-product 
of the low flow rates is then a low particle velocity in the 
bed and correspondingly a relatively low bed wall heat transfer 

coefficient (h ) . 

w 
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Ap (in. water ) A P <1n - water) A p (in. water) 



100r 



Mu 1 1 i -o r i f i ce 




Bed atr velocity (ft/sec) 




0.0625in. sheet metal with 
a 0.052in. hole per square 
inch of area. 




Bed air velocity (ft/sec) 



Screen 




300-mesh screen on 0.01?1n. 
thick plate, with twenty- 
four 0.021in. holes on 
triangular diamond pitch. 




Bed atr velocity <ft/sec) 



Po rous p 1 ate 




0.125in. s i ntered-metal 
plate of 0.0008 in. mean 
pore opening. 



Figure ?. Gas distributor effect on fluidization onset. 
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The general equation [36] for the minimum fluidization 



mass flow rate is 



J mf 



688 D p U82 { P f ( P s - P f ) }°* 94 
,, 0.8 8 



(3.3) 



with D , the average particle size, 
P 



D 



P 



1 




(3.4) 



Equation (3.3) can only be used, however, when the Reynolds 
number 

0.03445 D G c 
p mr 

Re m£ = (3-5) 

y 

is less than ten. Should the value of Re m £ be greater than 
ten, then the alternative is to use the nomograph [72] of 
Figure 8, with the calculated value of D p , y and Pf(p s “ Pf) 
to recalculate the value of G ^ . If, however, the value of 
D,p or p f (p “Pf) does not appear on the nomograph, then 
go to Figure 9 [72] and with the calculated value of Re m £, 
from equation (3.5), determine the correction factor to be 
applied to the value of G ^ calculated from equation (3.3) . 

The new value of G ^ is then calculated from the equation 

^ ^mf ^ new “ ^ ^mf ^correction ^ ^mf ^ calculated* (3.6) 
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NOTE-Left (or right) scale of (p s -pf)Of 
corresponds to left (or right) 
scale of G m ^ 
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Figure 8. Nomograph to calculate G m f. 
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Figure 9 . 



Correction Factor for Re>10. 
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It should be mentioned that there are numerous corre- 



lations proposed for the minimum fluidization mass velocity 
(Gmf) that are applicable to the experimental results of 
various researchers. These predictions are not in very close 
agreement and are a function of different apparatus construc- 
tions, fluids used, particles used, etc. 

One variation considered worthy of note is a result 
of the work accomplished by Wen and Yu [73]. They demonstrated 
that when particles of more than two different sizes are fluid- 
ized, they either separate completely, forming distinct layers 
or intermix forming one layer of rather uniformly distributed 
particles. For a bed containing two different particle sizes, 
if the particle diameter ratio is less than 1.3, the particles 
were observed to intermix. If, however, the particle size 
ratio was greater than 1.3, the smaller particles fluidized 
separately above the larger particles. When this occurs, 
the voidage must be obtained for each layer and special pro- 
cedures, delineated in Reference [73], used to compute the 
fluid mass velocity. 

5 . Bed Expansion 

This area of bed behavior is important in that it 
resolves the minimum bed chamber height that must be designed 
to enable fluidization to occur. Unfortunately, although 
gas-fluidization processes are more prevalent than liquid- 
fluidization processes, the mechanics of the gas-fluidization 
process is not as well understood in comparison to the liquid 
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process. It has been found that the liquid systems will ex- 
pand smoothly almost to a voidage of 1.0. With a gas system, 
however, somewhere along the course of the expansion, a point 
of mechanical instability is reached which degrades the expan- 
sion process. The location of this trip point or the reason- 
ing as to why it occurs has so far eluded researchers. There 
is some speculation, however, that it may be caused by the 
formation of aggregates of particles in the bed and that it 
is somehow related to the fluid velocity. 

Based on analytic and experimental data. Leva et al . 
[36,40] formulated a simple relationship between the fluid 
mass velocity and bed particle voidage that enabled the 
generation of a prediction model for bed expansion. Some 
examples of this gas-solid expansion data can be found in 
Figure 10. What is shown is that at higher flow regions the 
data will align along curves rather than linear lines and at 
intermediate flow regions the reverse is true. Note that as 
the voidage increases, (1 - e)/e 3 will decrease. Also, at 
elevated flow rates, the slopes (m) become less negative in- 
dicating that aggregation becomes more prevalent while the 
linear lines indicate little or no aggregation (intermediate 
flow region) . It is also noted that small-particle beds may 
be expanded over substantially wider flow ranges than larger- 
particle beds before aggregation conditions will set in. 

This is indicated by the fact that, with smaller particles, 

departure from the straight line occurs at higher values of 
the reduced mass velocity than with larger particles. 
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NOTE-Rgg reg at i on sets in where 
the data begins to deviate 
from the dashed lines. 
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Figure 10 



Gas-particle expansion data 



As far as practical methods of correlation are con- 
cerned, the data is most cases may be represented by straight 
lines. This will produce the relationship between slopes (m) 

and particle size (D ) indicated in Figure 11 [36,40]. This 

P 

correlation is based on virtually all of the available bed- 

expansion literature data to date. 

To calculate a predicted expanded-bed height of the 

gas-fluidized particles, two basic methods may be used. The 

choice of method will depend on the bed-particle diameter. 

However, a maximum mean particle diameter (D ) of about 0.015 

P 

in is recommended, otherwise the fluidization efficiency will 
be so low as to promote slugging. Fluidization efficiency 
(r|) is defined by the equation 

G f - G e 

0 = , (3.7) 

G f 

where and are respectively defined as the fluidizing 
gas mass velocity flowing through a bed at a given instant in 
time and the gas mass velocity required to initially expand 
the bed. Once slugging sets in the bed expansion height 
will begin to fluctuate as the gas bubbles expand and collapse. 
When this observed condition occurs no correlation to date 
has been found to give a predictable bed height, 
a. Method One 

(1) With D known, m may be found from Figure 11. 

P 

(2) Determine from Figure 3 or equation (3.2) . 
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Figure 11. Gas— f 1 u i d i zed slope values vs. particle diameter. 
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(3) Calculate G m ^ using equations (3.3) and/or 
(3.6) or Figure 8. 

(4) Draw fluidization line on log-log paper with 
slope m, see Figure 12 [40]. 

(5) For the chosen value of G^, find the corre- 
sponding value of £ ^ . 

(6) Solve the voidage equation, equation (3.2), 

for Lf by substituting and for and e^. 

b. Method Two 

(1) Calculate G^ using equations (3.3) and/or 
(3.6) or Figure 8 . 

(2) Calculate the reduced mass velocity, Gf/ G m f 

(3) Go to Figure 13 [36] with the reduced mass 
velocity and D to find the fluidization efficiency (n) • 

ir 

(4) Go to Figure 14 [36] with n and to evaluate 
the bed expansion ratio, R. 

(5) Calculate the static bed height L ^ from 
equation (3.2) . 

(6) Multiply L ^ by R to get the expanded bed 

height, . 

Once the expanded bed height has been calculated 
there still remains one check to be performed on its validity. 
This is a check on the extent of height fluctuations that may 
occur for the flow rate chosen. It has been found that under 
conditions of high gas flow rates that an oscillation of the 
expanded bed height can occur. This concept has been termed 
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Figure 12. Gas Fluidization and expansion lines. 
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Figure 13. Fluidization efficiency vs. reduced mass velocity. 
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Figure 14. Fluidization efficiency vs. expansion ratio. 
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the fluctuation ratio, r, the ratio of the highest to lowest 
levels which the bed top occupies for a particular gas flow 
rate. The value of r can be determined from the equation 

( G f " G mf } 

r = exp (3.8) 

G 

mf 

where slope m" is related to the particle diameter as shown 
in Figure 15 and is based on data of Lewis et al . [42] . The 

closer the fluctuation ratio to 1.0 the more accurate the 
expanded bed height estimate. 

B. HEAT TRANSFER PHENOMENON MODELS 

There is little agreement in the reported experimental 
values of the heat transfer coefficients in the literature. 

The reasons can best be summarized as follows: 1) a variety 

of contacting patterns between particles and the fluid re- 
sulted in different models from which to base the average 
heat transfer coefficient; 2) the temperature measurements 
taken in the beds by thermocouples were interpreted differ- 
ently (usually the readings were related to the gas temper- 
ature in some manner) ; or 3) that direct measurement of 
solid particle temperatures are unreliable with present 
state-of-the-art methods, especially for the smaller size 
particles (less than 100 microns) . 

As a result, due to the generally low accuracy of in-bed 
temperature measurements and the oversimplifications in the 
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Figure 15. Slope m' In relation to particle diameter. 
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flow models, heat transfer coefficients so far reported have 
varied widely. Although there is no complete theory that 
can predict the heat transfer data for fluidized beds, there 
are several theories that correlate the published data. Each 
resulting heat transfer model has some experimental or design 
situation for which it is appropriate. What follows is an 
overview of these various models all of which are based on 
cylindrical bed geometries. 

1 . Steady State Conduction Across the Gas Film 

In this model heat is conducted through the gas 
boundary layer near the heat exchange surface. In order to 
explain the high values of the heat transfer rates obtained 
experimentally, the gas layer is assumed to be scoured by 
solids moving along the heat exchange surface, thus decreas- 
ing the boundary layer thickness near the wall. This model 
was developed by Leva et al . [37], Dow and Jakob [13], and 

Levenspiel and Walton [41] . 

Dow and Jakob [13] studied the gas-fluidization mix- 
tures and classified them as either a dense or lean phase 
mixture based on the gas velocity. The dense phase occurred 
at gas velocities close to u m ^ while the lean phase was close 
to u t . Their model applied to the mechanically smooth dense 
phase fluidization although they did observe bubble formation 
and slugging in the dense phase region of fluidization. 

The model's motion. Figure 16, is primarily upward 
at the center of the bed and downward along the wall . 



56 



Temperature measurements showed the entire bed was at a con- 
stant temperature except for a small region at the bottom 
and a narrow layer near the wall. Heat is tranferred radially 
through a thin air film to the particles and air moving down- 
ward outside this film. The particles carry the heat to the 
bottom of the bed where thermal equilibrium is attained almost 
instantaneously with the incoming cold air. The bombardment 
of small particles prevents the formation of a laminar bound- 
ary layer, and leaves a very thin laminar sublayer and a 
thicker turbulent outer layer. 

The heat transfer rate from the surface to the fluid- 
ized bed was found to be independent of the thermal conduc- 
tivity of the solid particles, and the large isothermal 
region of the fluidized bed could be destroyed if the down- 
ward flow of the small particles was short-circuited into 
the main flow before they reached the bottom, i.e. under 
slugging conditions . 

Levenspiel and Walton [41] presented a similar model 
where the resistance to the heat flow is due to the laminar 
gas layer which is destroyed by solid particles passing 
through it. Thus the average thickness of the laminar layer 
is much less than in an empty tube. They assumed that the 
particles are of uniform diameter (d ) , stationary, and 
arranged in equally spaced horizontal layers. The boundary 
layer formation at the wall is similar to that of a flat 
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plate, see Figure 17. At points of contact with the station- 
ary particles, the boundary layer is destroyed and starts 
once again. 

Unfortunately, the two models do not take into account 
the influence of the solid particles on the heat transferred, 
so they must be considered of limited application. 

2 . Unsteady Heat Conduction by Single Particles in 

Direct Contact with the Heat Exchange Surface 

With the advent of photographic studies of the bound- 
ary wall during fluidization, it was shown that the particles 
tended to associate in groups which moved slowly along the 
wall or remained in contact with the wall a certain length of 
time. Those observations led Botterill et al . [6] to propose 

that heat was transferred from the vertical surface to the 
adjacent layer of particles arriving together at the wall. 

If the residence time of the particle near the wall was long, 
heat would penetrate further into the bulk of the particles 
and if the residence time was short, heat exchange calcula- 
tions would be limited to only one particle, the nearest to 
the surface. They confined their theoretical considerations 
to an isolated particle lying in contact with the wall. 

After much experimentation this model was extended 
to cover different particle sizes and materials and a wide 
range of residence times for the particles at the surface. 

The results showed the validity of the model for short resi- 
dence times with the observed effect that the heat transfer 
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Figure IS. Heat transfer mechanism and particle motion. 




Figure 17. Laminar boundary layer thickness between 
two 1 ayers . 
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coefficient is a decreasing function of the residence time. 
Again, this model can only be applied to mechanically smooth 
dense phase fluidization as the model's validity breaks down 
in the lean phase or when bubbling/slugging occurs. 

3 . Unsteady Transfer of Heat to "Packets" of Particles 
Which are Renewed by Violent Disturbances in the 
Core of the Fluidized Bed 

In an attempt to present a mechanism for the heat 
transfer in bubbling, dense phase beds, Mickley and Fairbanks 
[47] proposed a model of a small group of particles moving as 
individual packets. They based their model on the fact that 
the voidage of the dense phase bubbling beds was found to 
be close to the voidage of quiescent beds. 

In the dense phase fluidized bed, each particle was 
expected to be in contact with several neighbors most of the 
time. The packets are not considered permanent; they have a 
finite life expectancy; and their voidage, density, heat con- 
ductivity, and heat capacity are assumed to be the same as 
those of the quiescent bed. See Figure 18. 

This model assumes that as a packet at bed tempera- 
ture T^ comes into contact with a flat surface at a tempera- 
ture T that unsteady heat transfer will start upon contact, 
w 

Further, if A is the contact area between the nacket and the 
m 

wall, and the packet is homogenous, then the heat transferred 
after some time t will be given by the equation 



a . = h . T A 
-wi wiL m 



( T 



w 



- T b > 



(3.9) 
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where the local instantaneous heat transfer coefficient is 



given by the equation 



h 

wiL 



k p c 
m m m 



0.5 



TT t 



(3.10) 



Mickley and Fairbanks [47] also presented two devia- 
tions from their original model which covered the phenomenons 
of slug flow and side mixing. The first, on slug flow, was 
for fluidized beds with very low gas flow rates and no turbu- 
lent mixing. In this model, Figure 19, the particles were 
to move downward along the outside wall of the bed. The flow 
is not smooth but is spasmodic owing to the intermittent 
passage of ascending gas bubbles near by in the core of the 
bed. However, assuming that the particles moved downward at 
a constant velocity (u s ) , the age of all the packets at a 
distance l from the bed top was t = H/u and the resulting 
average heat transfer coefficient for a heater of length L„ 

li 

could be found from the equation 

0.5 

(3.11) 

The second model deviation, on side mixing, was for 
beds with large heat transfer surfaces that were highly 
turbulent. To represent this case, Figure 20, they assumed 
that as the particle packet is moving downward some of it 
is exchanged with some other packet at temperature T^, brought 



h = 2 
w 



k p c u 
m m m s 



7T L 



H 
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Figure 18. Transfer of packet to mail. 




Figure 19. Slug flow of particles over a surface. 
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in sideways from the core of the bed. The average number of 
replacement packets at the wall per unit time resulting from 
side mixing was then defined as s, resulting in the local heat 
transfer coefficient equation 



wL 



= ( k 



m 



m m 



s ) 



0.5 



(3.12) 



It should be noted that for these penetration type 
models the residence time of the particle is of fundamental 
importance and has generated a considerable amount of research 
both in its determination and its effect on the various models 
of this section. 

Patel [53], in a further spin-off of Mickley and 
Fairbank ' s original models presented two surface renewal 
models that can be characterized as follows. 

In model one. Figure 21, a packet initially at the 
bed temperature T.^ arrives at the wall (temperature T^) at 
time t = 0. The packet is assumed to have the same proper- 
ties as the bed at minimum fluidization. It receives heat 
from the wall through a contact resistance and returns to 
the bed after some length of time. During that residence 
time, heat penetrates a distance x; beyond this distance 
the temperature is constant and equal to the bed tempera- 
ture T^- 

In model two. Figure 22, the fluidized particles 
are assumed to be spherical and of the same diameter (d^) . 

A particle at the bed temperature T^ arrives at the heating 
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Figure 20. Down-flow of particles with side mixing. 
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Figure 21. Patel Model I. 
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surface (temperature T ) and receives heat from the fluid 
adjacent to the wall. This fluid is assumed to be stagnant 
and at wall temperature T . While the particle is at the 
surface it also loses heat by conduction to a packet of par- 
ticles of thickness x situated between the wall particle 
and the bulk of the bed. The packet is assumed to have the 
same properties as the bed at minimum fluidization conditions. 

Patel by no means had the only variation on Mickley 
and Fairbank's packet theory. Numerous other variations 
exist with different h^ values being correlated by each 
researcher. Appendices A through C give some of the result- 
ing correlations. 

4 . Steady Convection Through the Emulsion Layer 

van Heerden et al. [66] observed that the heat capa- 
city of the particles per unit volume is of an order of mag- 
nitude of thousands times greater than that of the gas and 
that the mean particle velocity is much lower than the gas 
velocity. Therefore, the largest portion of the heat will be 
transferred by the moving particles. They assumed that the 
gas temperature was that of the particles and as a result 
that the gas only provided a suspension for the particles. 

They observed negligible radial motion of the bed 
particles and therefore concluded that the radial heat trans- 
fer was a function of the thermal conductivity of the sus- 
pension. Experimentally, they verified that the heat transfer 
coefficient is large for short heat exchange surfaces and vice 
versa. 
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Wicke and Fetting [75] furthered van Heerden's work 
and proposed a model, Figure 23, for those observations. In 
their model, heat (<3 W ) from an exchange surface was first 
transferred by conduction through a gas layer whose thickness 
was 6 g . This heat was then divided into two components: 
namely heat (q z ) taken by the particles flowing parallel to 
the surface in a second zone of emulsion of thickness 6 , and 
heat (q r ) transferred into the core portion of the bed by 
interchange of the particles. In their calculations, how- 
ever, they neglected the term q r . 

The last model to be presented is actually a combina- 
tion of several heat transfer mechanisms proposed in earlier 
models . 

5 . Kunii and Levenspiel Model [34] 

In this model, Figure 24, four mechanisms operate 
together . 

a. A film of gas coats the surface. Its thickness 
(6 g ) is large or small depending on whether a bubble of air 
is near the surface or the emulsion is uniform and close to 
the surface. 

b. Some solid particles are in direct contact with 
the heat exchange surface. 

c. There is a layer of emulsion with thickness 5^ 
which flows along the wall. 

d. Part of the emulsion layer is replaced occasion- 
ally by fresh emulsion coming from the core of the bed or by 
bubbles rising along the wall. 
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Figure 22. Patel Model II. 




Figure 23. Wicke and Fetting Model. 
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Figure 24. Kunii and Levenspiel Model. 
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Depending upon the fluidization conditions and the 
position of the heat exchange surface in the bed, emphasis 
can be placed upon any combination of the above four 
mechanisms . 

C. FLUIDIZED BED VARIABLES THAT AFFECT HEAT TRANSFER 

In addition to the large number of models that have been 
proposed by various researchers to explain their experimental 
results in heat transfer, there are a number of variables 
which have been found to affect the heat transfer rate. In 
general, the factors which will control the heat transfer 
characteristics of the bed will be those variables which 
determine the fluidizing properties of the bed. In this 
section an attempt will be made to list some of the more 
important factors and their effect on heat transfer. An 
excellent indepth review of general fluidized bed character- 
istics was researched by Gutfinger and Abauf [28] . Addition- 
ally, Appendices A through C present a list of nondimensional 
correlations reported by different investigators while trying 
to correlate their experimental findings. 

1 . Properties of the Fluid 

a. Density (p^) and Specific Heat (c^.) of the Fluid 
The dependence of the heat transfer coefficient 
on the specific heat of the gas has not been clearly deter- 
mined as the existing data shows no common trend. However, 
Gelperin and Einstein [18] have predicted that there would 
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be an increase in the heat transfer rate with an increase in 



c^ or the product p f c pf when the fluid was at high pressures 
and velocities. 

b. Gas Viscosity (y) 

The majority of researchers agree that with an 
increase in the gas viscosity the wall heat transfer coeffi- 
cient will decrease. 

c. Thermal Conductivity of the Gas (k^) 

This is the fluid property that has the most pro- 
nounced effect on the heat transfer coefficient. Depending 

on which researcher is consulted the value of h will increase 

w 

proportional to k^ to the power of 0.5 - 0.67. 

2 . Properties of the Bed Particles 

a. Diameter (D ) and Shape of the Particle (<j) ) 

P s 

A significant number of researchers [37,47] have 
concluded that the heat transfer coefficient decreases when 
large diameter particles are used. However, the results of 
Sarkit et al . [60] as reported by Gelperin and Einstein [18] 
show that for the laminar flow region the heat transfer 
coefficient varies inversely with particle diameter and for 
the turbulent flow region it has a direct variation. 

The values for the heat transfer coefficient tend 
to increase as the sphericity (4> s ) increases (rounder and 
smoother particles) . 

b. Density of the Particle (p s ) 

Experimentation has shown that the wall heat 

transfer coefficient increases with the increasing density 
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of the solid particles. However, this increase slows in 
its magnitude as the fluid flow becomes more turbulent. 

c. Specific Heat of the Particle (c ) 

ps 

The heat transfer coefficient increases with an 
increase in the specific heat of the particle. 

d. Thermal Conductivity of the Particle (k s ) 

Based on measured experimental results over a 

wide range of particle thermal conductivities it has been 
determined that the effect of the thermal conductivity of 
the particle on the heat transfer coefficient is slight, so 
much so that it is usually neglected as a factor of concern. 

3 . Fluidization Conditions 
a. Gas Mass Velocity (G) 

The increase in the heat transfer coefficient is 
proportional to G n where n will vary according to the experi- 
ment being conducted. Its range has been shown to be quite 
large. About the only facts that have been agreed upon by 
researchers for this variable is that of the general shape 
of the gas mass velocity versus heat transfer coefficient 
curve [72], Figure 25, and that there is some maximum value 
for h w which exists at some gas mass velocity G £. Tables 
listing the nondimens ional correlations for both the rising 
and falling sides of the curve, for the maximum heat transfer 
coefficients, and for the special gas flow rates when the 
hwmax i s attained are presented by Einstein and Gelperin 
[14]. 
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Figure 25. Heat transfer coefficient vs. mass velocity. 
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b. Bed Porosity or Voidage (e) 

Based. on film theory the heat transfer coeffi- 
cient will vary inversely with some function of the bed 
porosity. That is, as the particle aggregate at the wall 
decreases so will the heat transfer rate. However, there 
has been some experimental data which has demonstrated the 
opposite effect. 

4 . Geometric Variables 

a. Bed Diameter (D) 

It has not been demonstrated that there exists 
a dependency of the heat transfer coefficient on the bed 
diameter . 

b. Diameter of an Immersed Heater (D u ) 

rl 

Researchers are in general agreement that the 
wall heat transfer coefficient increases with a decrease in 
the immersed heater (tube) diameter, and that the heat trans- 
fer rate diminishes as the diameter of the heater is increased. 
These effects are valid up to a particle diameter of 10 mm 
(0.39 in) after which h w is independent of the heater 
diameter . 

c. Fluidized Bed Height (L) 

There has been no dependence shown between the 
heat transfer coefficient and the fluidized bed height. 

d. Length of Heat Exchange Surface (L^) 

This is a variable that has had several differ- 
ent results on its effect on the heat transfer coefficient 
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all of which have been substantiated by additional research. 
Basically, one research group has found that the heat trans- 
fer coefficient decreases with an increase in the length of 
the heat exchange surface. Another research group concluded 
that if l h /D was greater than 7 for an external heat transfer 
surface, then the heat transfer coefficient was independent 
of heater length and that internal heat exchange surfaces had 
no observed effect on the heat transfer coefficient. And yet 
a third research group found that if the heat exchange surface 
was located above the active section of the bed (mixing 
region) , then the heat transfer coefficient would be inde- 
pendent of the heater's length. 

e. Heat Exchange Surface Roughness 

It has been determined [31,68] that the heat trans- 
fer coefficient is sensitive not only to the type and orient- 
ation of the roughness but also to the size of the particles 
relative to the size of the roughness. When the particles 
are of a large size relative to the roughness, the gas con- 
duction paths will be longer and the heat transfer will be 
lower as compared to a smooth tube or surface. See Figure 
26a and 26b. 

When particle and roughness sizes are similar, 
the multiple contact should reduce the mean gas conduction 
paths and gives rise to a considerable increase in the heat 
transfer. Figure 26c. However, possibly due to particle 
capture and the consequent increase in particle residence 
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Figure 2G. Contact geometry of surf ace-particle. 



75 



times, the increase in the maximum heat transfer coefficient 
for vertical grooves, is of the order of 50% for the smaller 
particles and 18% for the larger particles. The particle 
capture and inability of the bubbles to replenish particles 
in the case of horizontal grooves undoubtedly contributes to 
this difference. 

Where the size of the particles is much smaller 
relative to the roughness size, particle capture is not a 
serious problem, but there is some improvement in particle- 
surface contact as well as a larger effective area for heat 
transfer, Figure 26d. 

5 . Miscellaneous Variables 

a. Entrance Effects 

Entrance effects in a flow situation cannot be 
ignored or avoided. Control or limiting entrance effects can 
only be accomplished by careful design of the gas distributor, 
which in turn, will influence the configuration of the active 
(mixing) region of the fluidized bed. 

b. Temperature and Pressure of the Fluidized Bed 

The heat transfer coefficient is independent of 

the pressure of the fluidized bed or fluidizing fluid. The 
temperature of the fluid, however, is another matter. In 
general it can be assumed that the heat transfer coefficient 
is independent of the temperature of the bed or fluid. 

Studies by Zabrodsky [79] and Yoshida et al . [77], however, 

indicate that the radiation component of the overall heat 
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transfer coefficient must be taken into account when the bed 



temperature exceeds 900-1300 K. In turn, the magnitude of 
the radiation component effect will be dependent on the bed 
container geometry . 

D. IMMERSED BODIES 

Experimental work conducted in fluidized beds can be sub- 
divided into two distinct categories. The first deals with 
all experimental work that uses an external heat transfer 
surface, i.e. the wall of the fluidized bed container itself 
is used as the heat exchange surface. The second category 
uses an internal heat transfer surface. In this case, bodies 
of different shapes (sphere, cylinder, plate, coil, etc.) 
which are immersed in the bulk of the fluidized bed become 
the heat transfer surface. 

While the previously discussed variables and character- 
istic bed behavior are equally germane to either category, 
the use of immersed heating surfaces will have additional 
factors that will need to be taken into account during a 
design process. The location of the heat source with respect 
to the gas distributor plate and its plane of orientation 
will have a significant effect on the heat transfer coeffi- 
cient. The results of the variations are sometimes contra- 
dictory and depend strongly on the local hydrodynamic 
conditions in each case. 

Limited experimentation has been conducted in the use of 
spheres, flat plates and small cylinders with the data 
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obtained being sketchy at best. The bulk of the immersed 
body experimentation has been with tubes for which additional 
variations in heat transfer tube design (size, shape, spacing, 
gap, pitch and material) will complicate the general heat 
transfer mechanism. An excellent indepth coverage of this 
topic was conducted by Saxena et al. [61] and is summarized 
as follows . 

1 . Horizontal Tubes 

The total heat transfer coefficient (h ) for single 

w 

cylindrical tubes increases with factors that increase the 

gas thermal conductivity near the tube wall (i.e. temperature 

at the tube surface) and reduces the particle residence time 

(i.e. decrease in tube diameter and increase in fluidization 

velocity). However, at higher gas velocities the decrease 

in bed density reduces the heat absorption capacity of the 

fluidized particles, resulting in a decrease in the heat 

transfer coefficient (h ) . The heat transfer coefficient 

w 

also increases with decreasing particle diameter (dp) • 

The heat transfer coefficient is highest for tubes 
with a lenticular cross section for dense moving packed beds. 
Gas pockets are not created in the flow wake (downstream 
side) of a lenticular tube, and the particles are not stag- 
nant on the upstream side of the tube. 

The presence of surrounding tubes in an array affects 
the total heat transfer coefficient for an individual tube 
by resisting particle movement and by altering the particle 
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temperature from that of the bed bulk, which reduces the tem- 
perature driving force. This is at least one of the reasons 
why the heat transfer coefficients for the lower tubes are 
greater than for the upper tubes, and a close spacing of the 
tubes results in a decrease of the heat transfer coefficient. 

Fins will increase the rate of heat transfer to the 
tubes. However, their effectiveness increases, reaches a 
maximum, and then decreases with decreasing spacing distance 
and increasing fin length [54] . 

An increase in tube-to-tube spacing in a bundle of 
short-finned tubes, which reduces the resistance to particle 
movement between the tubes, increases the total heat transfer 
coefficient. However, an increase in tube-to-tube spacing 
for a bundle of long-finned tubes does not noticeably in- 
crease the heat transfer rate, implying that the resistance 
to particle movement between fins predominates. 

2 . Slanted Tubes 

There does not appear to be any particular advantage 
in slanting the heat transfer tubes. Minimum heat transfer 
coefficients have been observed at 45° for bare tubes and at 
60° for finned tubes with respect to the horizontal bed axis. 
The heat transfer coefficient values obtained were in compari- 
son to the values obtained for tubes of the same configura- 
tion but in a horizontal and vertical position [23] . 

3 . Vertical Tubes 

The total heat transfer coefficient for a single ver- 
tical tube generally increases with a decrease in tube length 
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and diameter, a decrease in the fluidized particle size (sizes 
less than about 1 mm (0.039 in) at atmospheric pressure) , and 
an increase in fluidization gas velocity up to u fc . 

The heat transfer coefficient for a tube within a 
bundle increases with an increase in the tube spacing and 
is greater for tubes located nearer to the center of the 
bundle, at least for small bundles. 

The heat transfer coefficient for a finned tube 
generally increases with the fin height and a decrease in 
the frequency of fins along the tube length (in the range of 
21-10 fins/in) . Heat transfer also improves with roughened 
and grooved tubes . 
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IV. EXPERIMENTAL APPARATUS AND PROCEDURES 



This section will describe the details and procedures as 
related to the experimental apparatus . The methodology and 
experimental procedures will also be discussed. 

A. DESCRIPTION OF THE FLUIDIZATION APPARATUS 

The fluidization container for the experimental investi- 
gation was constructed of one half inch thick clear plexiglas 
(Grade GM) . The inside dimensions of the rectangular enclo- 
sure were variable as one wall was capable of movement within 
the boundaries of the container. This movement allowed the 
variation of the inside dimensions from six inches wide by 
six inches deep by eighteen inches high to twenty four inches 
wide by six inches deep by eighteen inches high above the 
distributor plate. All joints were glued together and wood 
screws added in mounting the distributor plate for additional 
weight carrying support. Mounted within each of the two side 
walls is a five inche wide by ten inch high strip heater 
that has the heating surface on the same plane as the inside 
plexiglas wall. Backing the strip heaters are two layers of 
insulation (silicone rubber sponge and plexiglas) which 
provide both rigidity and a mounting platform for the strip 
heater. Static pressure taps and bed thermocouple taps were 
installed at evenly spaced intervals as shown in Figure 27. 
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Rir Distributor Tube 



Pressure drops were measured using the static pressure taps 
with the output going to a Meriam Type W Model A324 sixty 
inch water manometer or a Model 33KA35 thirty inch water 
manometer via a pressure manifold. Air and bed temperatures 
within the experimental apparatus are measured with copper- 
constantan thermocouples as described in the EXPERIMENTAL 
HEAT TRANSFER PROCEDURE AND APPARATUS section of this thesis 
and recorded on a Newport Model 267B-TF2, Multipoint (12 
Channel) Digital Temperature Recorder. The complete equip- 
ment listing for this experimental study is given in Appendix 
E. 

The distributor plate utilized in the fluidization con- 
tainer was constructed of one half inch thick plexiglas 
(Grade GM) with a grid of 0.125 in diameter holes bored at 
0.25 in center to center spacing. To cause the required 
pressure drop which aids the air distribution and to prevent 
the particles from falling into the air box, one layer of 140 
mesh stainless steel wire cloth was mounted on top of the 
distributor plate. 

The air was supplied by a Clements Cadillac Blower/ 
Suction Cleaner Model G12 with a nominal output of about 20 
scfm of air at a nominal line pressure of 1.5 psig. The flow 
output data is based on experimental testing of the blower 
since blower characteristic curves are not available. The 
air inlets to a two and one half inch inside diameter plexi- 
glas tube that has one half inch circular air ports spaced 
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at one inch center to center intervals along its length. The 
air ports are located at the underside of the tube allowing 
the air to exhaust downward into the air box. The plexiglas 
tube runs the length of the air box which has inside dimen- 
sions of twenty four inches wide by six inches deep by six 
inches high. The air is unfiltered and prevented from exit- 
ing outside of the variable bed inlet area by an air baffle 
mounted to the movable wall above the distributor plate. The 
schematic line diagram of the air supply system is shown in 
Figure 28. 

The air flow rate was measured with a Fischer and Porter 
Company Model 10A3565A Flowmeter (rated at 11.1 scfm) which 
was corrected for local barometric pressure . Photographs 
showing two views of the overall experimental apparatus can 
be found in Figure 29 and 30. 

B. EXPERIMENTAL HEAT TRANSFER PROCEDURE AND APPARATUS 

Four different rectangular bed configurations using a 
standard static bed height were investigated. For each bed 
configuration both an unheated and heated run were conducted 
with various pressure, temperature and visual data taken on 
each run. 

Once the bed size was adjusted, the bed of particles was 
initially fluidized at a point above minimum fluidization. 
When steady state conditions had been reached, the bed and 
air (inlet and outlet) temperatures were recorded and a 
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Figure 28. Diagram of Fluidization system, 
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Figure 29 . 



Photograph of experimental apparatus. 
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Figure 30. 



Photograph of experimental apparatus 
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comment recorded describing the visual condition of the exter- 
nal boundaries of the bed. This was the starting point for 
the run. 

As a separate study, bed pressure drop was recorded as a 
function of several flow rates . This information was used to 
determine the point of minimum fluidization for the particular 
bed configuration and to serve as an indicator for the pre- 
sence of slugging or channeling conditions within the bed. 

After the base data run was recorded, power was applied 
to the two Watlow Silicone Rubber strip heaters by way of a 
Lamba Regulated Power Supply, Model LK345A. The power supplied 
to the heaters was computed by measuring the voltage and 
amperage. Energy losses by conduction were accounted for in 
the heater configuration and an energy balance computed for 
each run. 

Since the purpose of this experiment did not include a 
need to determine a heat transfer coefficient, provisions 
were not made to determine heater wall temperatures or the 
temperature gradients close to the heater wall. On all runs 
the expanded bed height covered the strip heaters ensuring 
that all heat transfer was either to the bed or losses 
through the insulation. Radiation losses were neglected 
since bed/wall temperatures were considerably less than 1000 
K. 

At various heater power levels a series of temperature 
readings were recorded until flow characteristics were 
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determined for each bed configuration. The above experimental 
procedures were continued until all bed geometry configura- 
tions had been studied. 

The details of the heater configuration and mounting 
method are shown in Figure 31. Insulating end blocks were 
used for both mounting and end heat loss reduction for the 
two heaters. Heater insulation consisted of 0.5 in of sili- 
cone rubber sponge with a mounting block of 0.5 in of clear 
(Grade GM) plexiglas. Two insulated copper-constantan thermo- 
couples were applied to the outside of the plexiglas mounting 
blocks to measure their outer surface temperature. These 
temperatures were used to estimate the heat loss through the 
back of the heater by conduction for use in the heat balance 
of the model . 

The strip heater resistance element is constructed of a 
fine nicket alloy wire wrapped around a glass cord. The 
element wires are arranged in a parallel rectangular configur- 
ation thus confining any failure to a single element rather 
than the entire unit. Each heater is rated at 250 watts, 

115 volts with a temperature range of 0-200° F. A consider- 
ably higher temperature range is possible by cementing the 
heater to various metal plates, but is limited in this applica- 
tion due to the melting temperature of the plexiglas. 

To measure the bed temperatures a grid pattern of thirty 
retractable Omega Engineering ANSI T (Copper-constantan) 
twenty-four inch stainless steel sheathed exposed junction 
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Figure 31. Heater Rssembly. 
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thermocouple probes was used. By varying the depth of the 
individual thermocouples the boundary positions of the par- 
ticle and fluid flow cells could be determined and changes 
to the original pattern noted during the various runs. 

C. EXPERIMENTAL APPARATUS CALIBRATION 

Manometers used for static pressure drop and differen- 
tial pressure measurements were continuously rechecked for 
zero readings and readjusted as required. 

The temperature recorder and thermocouples were calibra- 
ted with a Rosemount Commutating Bridge Model 920A and 
Rosemount Calibration Bath Model 913A using liquid nitrogen 
and water. A least squares linear curve fit was applied to 
the resulting data for each thermocouple to obtain suitable 
constants for temperature correction. 

Technical data on the bed particles and experimental 
verification of the manufacturer's specifications can be 
found in Appendix F. 
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V. PRESENTATION AND DISCUSSION OF RESULTS 



Based on bed observations during the fluidizing process 
a very definite electrostatic charge built up within the bed 
container, see Figure 22. This phenomenon was characterized 
by the bed particles adhering to the walls of the bed after 
the air had been turned off and to the thermocouple probes 
during fluidization. Based on evidence provided by Miller 
and Logwinuk [49] this phenomenon will to some extent con- 
tribute to erratic results in the heat transfer and data 
measurements. Some of the factors which are thought to 
affect the magnitude of the electrostatic charge effects 
are (1) the dielectric constant of the bed material used, 

(2) the contact potential difference, (3) the ratio of fluid- 
ized bed particles to gas, (4) the properties of the gas, (5) 
the temperature of the solid-fluid system, and (6) the nature 
of the surface of the bed particles. For the purpose of this 
thesis, the electrostatic charge existance will be acknow- 
ledged but no data correction will be considered warranted. 
Further study on this area is, however, considered worthwhile. 

A. FLOW PATTERNS 

The particle movement along the wall at any level is in- 
fluenced by the movement within the fluidized bed at that 
level, as it is the movement of the inner particles which 
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drag the outer particles along the wall. Therefore the wall 
velocity profile should represent approximately the average 
particle velocity at the various levels within the bed. Had 
the motion been smooth with all particles starting at the 
top of the bed and flowing straight down to the bottom of the 
bed then the ideal, Figure 2, flow pattern would have occurred 
and the velocity of the particles could be assumed to be con- 
stant at the wall with a positive linear gradient into the 
axis of symmetry for the bed. In actuality this was neither 
observed nor anticipated within the rectangular bed under 
study due to the geometry involved. 

In the initial geometric configuration, that of a square 
column (6x6x12.5 in), the particle motion appeared to conform 
to the particle circulation pattern of Figure 33. The fluid- 
ization somewhat resembled waves in that there were intermit- 
tent fronts moving smoothly up through the bed which terminated 
in a percolating (rupture) action at the top of the bed. The 
domed bulge, Figure 34, of particles was observed roughly in 
the center of the column and after rupturing, the bed particles 
would be thrown to the sides where they would start down the 
bed walls. 

It was not possible to observe any distinct points of 
emergence into the bed of the particles in the downward tran- 
sit but a definite entry point was visible at the bottom of 
the bed. The characteristic shape of the non-movement area was 
elliptical along the bed bottom and curving up at the corners, 
see Figure 35. 



93 




Figure 32. Electrostatic charge. 
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Figure 33. Proposed particle cell circulation. 
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Figure 34. Square fluidized column — top. 
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Initially, the shape of the non-movement area was not 
elliptical. If the side of the apparatus as shown in Figure 
27 is considered to be the front face then bubbles were ini- 
tializing in the front left corner of the bed along the 
movable wall. This shift in symmetry of the fluidizing cell 
generation was resulting in a shallower non-movement area at 
the left front corner and a larger build up at the right 
front corner. The cause of the problem turned out to be 
that the distributor plate was not at a right angle to the 
air flow. Downward motion was observed at all four walls and 
the non-movement zone was located within one inch of the 
bottom of the bed. 

It was readily observable that the velocity of the par- 
ticles emerging at the center of the bed was higher than that 
transiting down the face of the wall. This differential 
while not empirically measured during this study was attri- 
buted to frictional drag at the wall both mechanical and as 
a result of the electrostatic charge present at the bed wall . 
There is some justification provided by Toomey and Johnstone 
[62]. Marsheck [43] also gives an excellent treatise on the 
velocity profile across a gas-fluidized bed. 

The height of the domed bulge would pulse during the 
fluidization more so as the air velocity was increased. If 
the air blower was fully energized with the initial bed con- 
figuration then the entire top of the bed would rise up as 
a mass and then collapse on itself. This was an example of 
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slug flow where the rising bubbles had coalesced into one large 
bubble stretching from one wall to another. Slug flow was not 
achieved or expected in any of the other configurations as the 
bed was not deep enough nor was the air velocity high enough 
to provide the right conditions for the bed cross-sectional 
areas present. 

For the air flow that just permitted particle circulation 
the bed height fluctuations were negligible although the par- 
ticles at the wall still moved downward in an irregular 
slightly surging movement. It is felt that at some point 
in the bed for this minimum fluidization velocity that the 
flow pattern internally is such that the flow pattern dis- 
tributes the particle velocities more evenly throughout the 
column. In other words, the circulating motion caused by 
the higher velocity in the center and slower velocity along 
the wall gives way to a more random or eddying motion that 
evenly distributes the small air bubbles within the two 
velocity boundaries thus preventing rapid bubble coalescence. 
These eddy cells once formed in the mixing region of the bed 
will be moved up through the bed by the surges resulting 
from the cell rupture at the bed top. Between each surge 
line it is theorized that only one eddy cell will exist and 
that it will grow as it progresses through the bed until it 
coalesces with the gas in the central cell core at the bed 
top prior to cell rupture. 
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It was also noted that as the thermocouple probes were 
lowered into the bed that a flattening of the bed top would 
occur. This effect became more pronounced as the probe pene- 
tration increased and is felt to be as a result of the probe 
providing a grounding channel for bubble/particle flow. In 
effect the probes acted as a calming agent in the agitation 
of the bed without changing the overall bed flow cell 
configuration . 

When the movable wall was positioned for a 6.25 inch width 

(A =37.5 in 2 ) and also for an 8 inch width (A = 48 in 2 ) the 
c c 

bed expansion observed at the minimum fluidizing air velocity 
was 0.125-0.375 in. When the bed width was then increased 
to 10 inches (A c = 60 in 2 ) and subsequently 12 inches (A = 

72 in 2 ) the bed expansion observed did not exceed 0.25 in. 

This measurement does not take into account the height in- 
crease of the cell dome prior to rupture. There was no 
general horizontal movement of the particles near the wall 
except at the top of the fluidized bed and in a rather 
turbulent mixing region near the bottom of the fluidizing 
bed. 

As the movable wall was repositioned thus increasing the 
bed air inlet area the bed became an aggregative fluidized 
bed. This resulted in multiple cell formations within the 
bed vice the single cell observed in the initial bed config- 
uration. Unfortunately, there was no observable pattern 
to the cell formations as their overall positioning was 
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Figure 36. 
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72 in 2 ). 
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random in nature, see Figures 36 and 37. Only one cell was 
observed to break the surface of the bed at any given instant. 
Occasionally, small bubbles were noticed breaking the surface 
of the bed between the larger bubble cells. This was thought 
to be caused by channeling between the larger bubble cell 
boundaries . 

There was no noticeable change in the outward appearance 
of the particle cell flow patterns in any of the bed configura- 
tions when the external heaters were energized during 
fluidization . 

B. PRESSURE DISTRIBUTION 

Static pressure readings were taken at various intervals 
along the height of the fluidized bed on all experimental runs. 
In order to determine the gas mass flow rates for minimum 
fluidization, measurements of the pressure drop across the 
bed were also made at various flow rates. All data obtained 
during the various runs can be found in Appendix H . The mass 
flow rate at the maximum pressure drop was then taken as the 
point of minimum fluidization in each case. The maximum 
pressure drop was used because it is the most sharply defined 
point on the plots, see Figure 38. 

The results of plotting the pressure drop versus gas mass 
velocity profiles produced very similar patterns which are in 
general agreement with the expected curve for the distributor 
design used. As can be seen from Figure 38, when the air inlet 
bed cross-sectional area was increased the curve would shift 
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Figure 38. Bed pressure drop vs. gas mass velocity. 
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to the left. This displacement may have been as a result of 
the lessening in the contribution of the mechanical wall 
friction drag to the total overall drag experienced by the 
bed with the increased air flow for fluidization. Neverthe- 
less, it would appear as the bed air inlet cross-sectional 
area was increased the gas mass velocity required for fluid- 
ization would also increase. Additionally, in all cases ob- 
served, the log Ap value for minimum fluidization was 
approximately 1.25. 

The results of the pressure measurements seemed to indicate 
that mechanically smooth fluidization was taking place. 
Channeling and/or slugging was not occurring at the minimum 
fluidization operating conditions although it was possible 
to demonstrate both channeling and slugging for the initial 
bed configuration. It was noted that during fluidization 
small fluctuations were visible in the liquid levels in the 
manometer that was measuring the static pressure in the bed. 
This operating characteristic is in agreement with similar 
conditions observed by Parent et al . [52] in their study of 

gas fluidization. 

The chain of events in measuring the pressure drop is as 
follows. Once the bed was initially settled the air flow 
was initialized to the bed. As the air flowed through the 
bed, a corresponding pressure drop occurred. This pressure 
drop increased linearly with the gas flow until at some point 
the bed began to expand. At this point the pressure drop 



102 



became nonlinear yet still increased to some maximum value. 
This maximum value was picked as the minimum fluidization 
point as previously stated. The pressure drop would then 
decrease until the bed started to circulate or fluidize at 
which time the pressure drop would stabilize at some point 
lower than the maximum value recorded earlier. The pressure 
drop remained steady as the air velocity was increased further 
but soon started to rise again. This new rise in the pressure 
drop across the bed was an indication that slugging was start- 
ing to set in and was especially prevalent in the initial bed 
geometry. It was also noted that this secondary pressure 
drop increase was less prevalent as the bed cross-sectional 
area was increased. 

C. TEMPERATURE DISTRIBUTION 

Air entered the air box at temperatures between 70° and 
120 °F and in every case left the bed at a temperature within 
1°F of that of the top of the bed, as measured at a point in 
the axial center of the bed. Only two data runs were made 
with the bed instrumented to record the temperature distribu- 
tion. The first run (A c = 37.5 in 2 ) had an average air inlet 
temperature of 10 8° F ±2.5°F and attained a constant air outlet 
temperature of 125°F ±0.4°F, see Figure 39. The second tem- 
perature instrumented data run ( A = 72 in 2 ) had an average 
air inlet temperature of 116° F ±2°F and attained a constant 
air outlet temperature of 129. 5°F ±0.7°F, see Figures 40A and 
40B . 
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In the graphics of Figures 39 through 40B it was necessary 
to plot the profiles for each thermocouple probe as there were 
too many variables to produce a true three dimensional temper- 
ature profile. The x-axis of each graph is the height of the 
probe tip in inches above the bed distributor plate. The y- 
axis of each graph is the temperature in °F. A probe grid is 
provide for each bed configuration with the probe number for 
reference. In all cases the temperature recorded by each 
probe above the 1.5 inch height was within the constant air 
outlet temperature range specified above for each data run. 

All temperature data recorded can be found in Appendix H . 

As can be readily observed, the bed quickly attains a 
constant temperature except for a very small region at the 
bottom of the bed and a narrow layer near the container walls. 
The area near the walls is by assumption as the experimental 
apparatus was not instrumented to read that area. 

It appears that the particles were heated by the heater 
as they transited the heater surface. Then, when they reached 
the bottom of the bed, thermal equilibrium with the incoming 
cold air was attained almost instantaneously. Except for the 
layer along the walls of the bed and in the small turbulent 
mixing region at the bottom of the bed, thermal equilibrium 
between the particles and the gas very nearly exists through- 
out the fluidized bed. 

Since the inside walls of the bed container were not 
instrumented with thermocouples it was not possible to 
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Figure 40R. Bed temperature vs. bed height (fl c = 72 in 2 ) 
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determine if there existed a temperature differential between 
the walls which did not have heaters and the bed bulk 
temperature . 

In order to determine the effect of the air velocity on 
the temperature readings of the thermocouple probes several 
runs where conducted at different air inlet velocities with 
the probe tip placed six inches above the top of the distri- 
butor plate. Small temperature fluctuations of at most 1°F 
were observed while the bed was fluidized and being heated 
under steady state conditions. The air flow was then ter- 
minated while monitoring the bed temperature readings. It 
was observed that the thermocouple reading would initially 
stabilize at a temperature within the range of the previous 
fluctuations and did not change for several seconds. A 
gradual rise then would commence due to the heat conduction 
from the heaters in the walls. Since the temperature rise 
did not exceed 1°F at the end of four minutes it would seem 
to indicate that the bed thermocouple readings were basically 
independent of the fluid velocity. It would also tend to 
support the assumption made earlier in this thesis that the 
particle temperature and gas temperature at any point in the 
bed can be considered to be at equilibrium. 

The only change in the mixing cell pattern within the 
fluidizing particle cell appeared to be the production of 
the temperature gradient at the bottom of the bed and an 
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increased bulk bed temperature. This was in comparison to 
a steady state bulk bed temperature equal to the air inlet 
temperature prior to energizing the plate heaters. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

The following conclusions may be made from the results 
of this investigation: 

1. Upon fluidization of the bed an electrostatic charge 
will build up within the bed. Although no data correlation 
error is used in this thesis, there most probably is some 
effect on data results. 

2. With a square column bed the particle flow cell formed 
is basically as shown in Figure 33. Only one major cell is 
formed with its axis of ascent roughly located at the center 
of the bed. 

3 . With an expanded rectangular bed several particle 
flow cells form in a random fashion. There is no pattern 
to the overall cell coverage and the bed can be termed a 
bubbling bed. 

4 . That since the bulk bed temperature remains constant 
throughout the bed that there must exist eddy cells within 
the bed that continually mix the bed particles as they move 
up through the bed. Otherwise, a gradient would form between 
the heater walls and the center of the bed. 

5. That the velocity in the center of any cell was 
greater than the velocity of the particles being forced down 
the walls. In the expanded bed the particles moving down the 
walls were a mixture of the various ruptured cells. 
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6. That mechanically smooth fluidization was not taking 
place at all times. Occasional instances of channeling were 
present and could be distinguished by small pertubations 

to the bed surface that did not change position - a simmering 
appearance . 

7. That the distributor plate used did provide enough 

of a bed pressure drop to allow fluidization to occur but that 
the tendency for channeling or slug flow was high if the air 
velocity was increased much past the minimum fluidization 
velocity . 

8. That a temperature equilibrium between the incoming 
fluidizing air and the bed particles was reached almost in- 
stantaneously in a small mixing region of the bed located at 
the bottom of the column. 

9. That the only temperature gradient that existed in 
the bed was located at the bottom of the bed where the mixing 
region was located. The remainder of the fluidized bed stayed 
at a constant bulk bed temperature irregardless of the number 
of flow cells present. 

10. That the air velocity did not affect the bulk bed 
temperature readings. It does, however, affect the depth of 
the mixing region at the bottom of the bed. 

11. That energizing the plate heaters did not visibly 
alter the particle cell flow pattern but that it did produce 
a gradient in the temperature profile of the bed interior. 
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12. That the rectangular bed configuration under study 
would trip into fluidization when the log Ap across the bed 
approached 1.25. 

B. RECOMMENDATIONS FOR FURTHER STUDIES 

Based on the results of the experimentation to date, 
research conducted in preparing this thesis and the problems 
experienced in trying to get the fluidization apparatus to 
work the following areas of study are recommended for future 
thesis study: 

1. A determination of what effect the electrostatic 
charge which builds up during the fluidization process has 
on the heat transfer or flow cell formation. 

2 . A study of the resulting changes in the cell formation 
patterns when bed particles of sizes larger than 1 mm are used. 
Particles in this size range are of much more practical use in 
the metallurgical industry. 

3. A study of the flow cell patterns with beads smaller 
than 300 microns which should give enhanced mechanical 
fluidization . 

4. A study of the effect of particle stratification 
(elutriation) when particles of different sizes are combined 
in the bed. The combination can be of varying sizes/rough- 
ness/sphericity or any other physical property of the bed 
particles which would promote the separation of the particles 
into distinct layers during fluidization. 
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5. A review of the cell pattern formation when using 
particles from various manufacturers (same or different bead 
materials) while maintaining a constant particle size. 

6. Conduct a pressure study on the bed with different 
distributor configurations . This would include such designs 
as porous plates, baffles, calming chambers, filter screens 
or various combinations of the aforementioned designs. 

7. Study the effect of changing the location of the heat 
source from opposite wall external heaters to other configura 
tional combinations such as all four walls, adjacent walls, 
single wall, submerged heater, etc. 

8. Conduct a velocity profile study similar to that con- 
ducted in Reference [4 3] . That study used a modified thermis 
tor anemometer probe to determine the flow pattern within a 
circular column bed. 

Additionally, prior to conducting future research the 
following apparatus modifications are suggested: 

1. Replace the left side sponge rubber gasket on the 
movable wall with a new gasket . 

2. Implant a pair of thermocouples in the unheated walls 
of the fluidized bed. 

3. Mount a pair of thermocouples on the surface of each 
strip heater to enable the monitoring of the heater wall 
temperatures . 

4 . Manufacture a stand for the cabling to the thermo- 
couple probes. The stand should be of sufficient height 
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to prevent the bending of a probe from the weight of its 
cable when the probe is retracted from the bed. 

5. Change the mounting method of thermocouples 35 
through 38 so that they are mounted in the backing plate 
vice on the backing plate. 

6. Unless the area of study is to focus on the electro- 
static charge built up during fluidization, provide some 
means to ground out this phenomenon. 

7. Provide a leveling system for the fluidization appara- 
tus that will ensure that the fluid flow direction out of 

the distributor plate is parallel to the bed container walls. 

8. The probe taps on top of the fluidizing apparatus 
should be bored to 5/64 inch for ease in inserting the probes 
during data runs . Additional nylon bushings should also be 
ordered. 

9. A second flowmeter of greater capacity than 11.1 scfm 
should be installed in the air supply line. This will be re- 
quired to measure the higher air velocities needed to fluidize 
the bed with further expansion. 

10. A cooler air source should be used in future runs to 
enhance the temperature differential at the bottom of the bed 
during fluidization with the wall heaters energized. 

11. A new cap for the fill vent is requird to prevent the 
glass beads from getting between the cap and plexiglas thus 
jamming the cap. 
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C . CONCLUDING REMARKS 



Particle motion is a complex function of the particle and 
gas properties, fluidization velocity, gas distributor and 
geometry of the bed. While a simple correlation can easily 
be made it must be tempered in its general application by the 
fact that it will most likely not give adequate allowances 
for the way in which the particles circulate and subsequently 
change with changes in the bed particle material, fluidizing 
conditions, scale and design of the fluidization apparatus. 

The interpretations presented in this thesis are based 
on the premise of visible particle movement at a given bed 
height and must be qualified by the fact that they represent 
a very limited amount of data. 

There is a general opinion among researchers that the 
bubbles of air passing through a fluidized bed grow in size 
by coalescence of smaller bubbles into larger bubbles. Even- 
tually given a deep enough bed the bubbles should coalesce 
into one large bubble located in the center of the bed in 
order to maintain the uniformity of the bed [43] . While 
this opinion has been generated as a result of experimenta- 
tion using circular beds it would appear that it can be 
equally applied to rectangular beds. No doubt had the 
expanded rectangular bed been deep enough, eventually it also 
would have conformed to the aforementioned opinion. 

The determination of particle flow patterns in a fluidized 
bed would be of extreme value in the establishment of a kinetic 
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or mass transfer model of a fluidized bed. The fact that 
there are a multitude of interrelated variables involved, 
unfortunately, has made the modeling a theoretical tool vice 
the needed operational scale-up tool required for industry. 
More study will undoubtedly be required until a workable 
industrial model can be produced. 

It has also been an observation that while the majority 
of the fluidization research was conducted by American 
researchers in the the late forties and early fifties, that 
the bulk of the recent research has been conducted by the 
Russians with some notable exceptions. This trend is neither 
understood nor considered warranted by the researcher con- 
sidering the very pronounce heat transfer improvement demon- 
strated through the use of fluidization. 

% 



116 



APPENDIX A: NONDIMENS I ONAL CORRELATIONS 

FOR EXTERNAL WALLS 



These correlations are but a sampling of the many 
fluidization correlations for external walls that have 
been theorized. 

1. Investigator - Baerg et al. [3] 

Correlation - 



h = h - 55 exp ( -0.012 { G - 0.71 p , 

w wmax c sb 



) 



Remarks - See Appendix C for h 
wmax 

2. Investigator - Bartholomew and Katz [5] 
Correlation - 

h w 1.56 + In ( Re Ga - 0,6 6 7 - 0.0 12 ) 

c - G -0.227 Pr 0 ’ 6 6 7 Ga 0 * 42 

Pf 



Remarks - None. 

3. Investigator - Brazelton [7] 

Correlation - 

D • G ' - °- 8 7 

( 1 - e ) " 1,3 

y 

Remarks - None. 



w 



c - G 1 
Pf 



= 0.72 
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4. Investigator - Das and Sarkar [11] 



Correlation - 



h D 
w 



0.2 4 



= 140 r 0,18 Re 0,09 



Remarks - None. 

5. Investigator - Dow and Jakob [13] 
Correlation - 



h D 
w 

— n ^ 


D 


0.6 5 


D 


0.1 7 


'< 1 - e > p s c ps 


0.2 5 


D G 


k f 


L 




d 

P_ 




£ p f c pf 




y 



Remarks - None. 

6. Investigator - Gamson [17] 
Correlation - 



St Pr O ' 6 6 7 = 2.0 Re ~ °* 6 9 ( 1 - e ) " °- 3 0 

m 



Remarks - None. 

7. Investigator - Lemlich and Caldas [35] 
Correlation - 

Nu = 0.055 Re 

Remarks - Low velocity liquid-fluidization. 
Correlation - 



Remarks 



St Pr 0-6 6 7 



D 


0.7 9 


l 


d 




Re . e 


P_ 




t 



High velocity liquid-fluidization. 
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8. Investigator - Leva [36] 



Correlation - 



Nu = 0.5 



d G 
P 



Remarks - None. 

9. Investigator - Leva and Grummer [37] 
Correlation - 



h = 3000000 k , D 
w r p 



D P G f n 



u R 



0.6 0 



Remarks - None. 

10. Investigator - Levenspiel and Walton [41] 
Correlation - 

' ^ 0.5 0 



h D 
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w r 


0*3 
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. Y j 


i i 
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h D 
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= 10 
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h 0.6 ( d G ) 0,7 0 

w p 



C r; G 

pf 
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10. Remarks 



D 



r 



tt d 

P 



6 ( 1 - £ ) 



Y = ( { 1 + a 2 } 1,5 - a 3 ) 



a = 0.041 ( 1 - e 




0.5 0 



u 



G 

G e = - 

£ 

11. Investigator - Mickley and Trilling [48] 
Correlation - 

0.0118 p G 0 - 263 
m 



Remarks - None . 

12. Investigator - Rao and Kaparthi [56] 
Correlation - 



Nu 



0.014 Re 



1.1 



1 - £ 



0.4 0 



Remarks - Semi-fluidized bed. 
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13. Investigator - Richardson and Mitson [57] 



Correlation - 



Nu = 119 Pr 



0.4 0 



c 


0.2 8 


’k 


ps 




s 


. c p f . 




. kf . 



0.0 2 



Re 



N 



D G 



0.0 7 5 



Remarks - Liquid-fluidized system with 



N = 0.02 



3.45 + 



14. Investigator - Rukenshtein [58] 
Correlation - 



Nu = 0.326 Re 0 ' 1 * 23 Pr 0 - 33 Ga oa 4 



Remarks - None. 

15. Investigator - Toomey and Johnstone [62] 
Correlation - 



log 



h 



w 




0.575 log 




+ 0.130 



Nu 



3.75 




Remarks - None. 
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16. Investigator - van Heerden et al. [66] 



Correlation - 



Nu = 0.58 



"b d g' 
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0.4 5 


p mf C ps 


0.3 6 


’ p f ‘ 
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_ p f c pf 




p mf 



0.1 8 



Pr 



Remarks - None. 

17. Investigator - Wen and Leva [72] 
Correlation - 



Nu = 0.16 



c s P. V* 9“*' 


0.4 0 


1 

CD 


k f 




R 



0.3 6 



Remarks - None. 



0.5 0 
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APPENDIX B: NONDIMENS IONAL CORRELATIONS 

FOR IMMERSED BODIES 



These correlations are but a sampling of the many 
fluidization correlations for immersed bodies that have 
been theorized. 

1. Investigator - Andeen and Glickman [2] 
Correlation - 



h D m 
w T 



G D t y 



0.3 2 6 



=900 ( 1 - e ) 



p c d 3 p g 
f p s ^ 



Pr 



0.3 



Remarks - Sometimes called the modified Vreedenberg 
correlation. 

2. Investigator - Gamson [17] 

Correlation - 



St Pr °-6 6 7 = 2.5 2 Re " °- 8 0 ( 1 - e ) " °- 3 0 

m 



Remarks - Vertical tubes . 

3. Investigator - Gelperin and Einstein [18] 
Correlation - 



Nu = 0.7 5 Ga°* 22 



T 

P 



0.1 4 



Remarks - Vertical tube bundle. 
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4. Investigator - Gelperin et al. [20] 
Correlation - 



£ 


= 0.7 3 Re 0 - 32 


1 


6 ( 1 - £ ) 


6 ( 1 - £ ) 


- 




_ 



Nu 



Remarks - Horizontal tube. 



0.3 2 



5. Investigator - Glass and Harrison [24] 
Correlation - 



Nu = 57.4 



’d n 


0.9 1 




0.1 
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N °* 3 4 
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_ D T_ 




u c 3 D 
mf M s 


d t 



o.o s 



u 



Remarks - Staggered tube bundle- N = 

6. Investigator - Huntsinger [29] 
Correlation - 



u 



mf 



Nu = 1.4 Re"' 3 J Pr 



0.5 3 r>„- 1.6 



Remarks - Coil. 

7. Investigator - Kruglikov [33] 
Correlation - 

Nu = 0.66 Re 0 - 86 Pr 0,29 



Remarks - Coil. 
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8. Investigator - Miller and Logwinuk [49] 



Correlation - 



h 



w 



1. 5 



G 0 - 5 



k 



s 



0.0 7 2 



k 



f 



2.4 



d °- 96 
P 




u 



0.8 



Remarks - Vertical tube bundle. 

9. Investigator - Petrie et al. [54] 
Correlation - 



h w °T 



= 14 



mf 



0.3 3 



Pr 



0.3 3 



0.6 6 7 



Remarks - Staggered and in-line tube bundles. 

10. Investigator - Sarkits et al. [60] 

Correlation - 



Nu 



0.013 3 Re 0 * 4 Ga 0,2 7 Pr 0 ' 33 
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Remarks - Coil. 

11. Investigator - Traber et al. [63] 
Correlation - 



Nu 



0.024 Re 0,65 Ga 0 - 10 Pr 0 ' 3 3 




o.i 



6 



Remarks - Coil. 
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12. Investigator - Vreedenberg [69] 



13. 



Correlation - 



Nu = a 
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mf mf 



0.3 5 



Remarks - Tube position varies. 
Investigator - Vreedenberg [70] 
Correlation - 
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G D p ( 1 - e ) 
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Pr 



0.3 0 



Remarks - Horizontal tube with- 
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is greater than 2050 but 



less than 2550 the value of Nu will be 
the average of the two Nu correlations 
above. 
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14 . 



Nu = 



15. 



Inves tigator - Wender and Cooper [74] 
Correlation - 



0.0119 ( 1 - e ) 
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Remarks - Vertical tube. 
Investigator - Zenz and Othmer [80] 
Correlation - 
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Remarks - Horizontal tube. 
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APPENDIX C: NONDIMENS IONAL CORRELATIONS FOR THE 

MAXIMUM HEAT TRANSFER COEFFICIENT 



These correlations are but a sampling of the many 
fluidization correlations for maximum heat transfer coef- 
ficients that have been theorized. 

1. Investigator - Baerg et al. [3] 

Correlation - 

0.00037 p sb ' 
d 

P 



h trTY ,,^ = 49 log 
wmax 



Remarks - None . 

2. Investigator - Chechetkin [8] 
Correlation - 

0.0017 Re . 0,8 Pr 0,4 
opt 

Nu = 

max , o.6 9 

P 



Remarks - Re . = 0.209 Ga 0,52 . 
opt 

3. Investigator - Chekansky et al. [ 9 ] 
Correlation - 
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wmax 



28.2 p 0,2 k °- 6 d "°- 36 
r s f p 
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Remarks - Staggered tube bundle. 
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4. 



Investigator - Gelperin et al. [19] 



Correlation - 

“r(' 

Remarks - Staggered horizontal tube bundle. 
5. Investigator - Gelperin et al. [21] 
Correlation - 




Nu 



max 



= 0.74 Ga 



0.2 2 



Nu = 0.64 Ga 
max 



0.2 2 



P 

D n 



0.0 9 



Remarks - Vertical tube bundle. 

6. Investigator - Gelperin et al. [22] 
Correlation - 



Nu 



max 



0.79 Ga 0 - 22 



1 




0.2 5 



Remarks - In-line horizontal tube bundle. 
7. Investigator - Jacob and Osberg [30] 
Correlation - 



h =h ( 1 - e ) ( 1 - exp { -pk £ } ) 

wmax o ^ r f 

Remarks - h Q and p determined by experimentation. 
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8. Investigator - Maskayev and Baskakov [44] 



Correlation - 

Nu =0.21 Ga 0,32 
max 

Remarks - Vertical cylinder in coarse particle bed 
with 1.4x10 s < Ga < 3x10 8 . 

9. Investigator - Pillai [55] 

Correlation - 

Nu = 0.36 5 0 0,8 2 Ga 0 - 22 
max 

Remarks - Sphere with Bi < 0.04 
10. Investigator - Ruckenstein [58] 

Correlation - 



Nu = Re 8,423 Ga 0,14 Pr 0-33 



max opt 



Remarks - Re , = 0.09 Ga 0,58 
opt 

11. Investigator - Sarkits [ 60] 
Correlation - 



Nu = 0.0087 Ga 0,42 Pr 0,33 
max 



1 

o 

*0 

_ _ 1 


0.4 5 


D 


0.1 6 


’h 

o 


1 

o 

Mi 

L 




d 

P 




d 

P 



0.4 5 



Remarks - Re ^ = 0.2 Ga 0 ' 5 (laminar region). 
Correlation - 



Nu =0.019 Ga °‘ s Pr 0,33 
rnsx 



~c 

ps 


0.1 


D 


0.1 3 


H o 


1 

4-1 

Ch 

o 

1 




_ d P_ 




d 

P_ 



o.i 6 



Remarks - Re ^ = 0.66 Ga 0,5 (turbulent region) 
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12. Investigator - Traber et al. [63] 
Correlation - 

Nu = 0. 021 Ga 0 ’ 4 Pr 0,3 3 
max 

Remarks - Re , = 0.55 Ga 0,5 . 

opt 

13. Investigator - Varygin and Martyushin [67] 
Correlation - 

Nu =0.86 Ga 0 ’ 2 
max 

Remarks - Re ^ = 0.118 Ga 0,5 . 

opt 

14. Investigator - Zabrodsky [78] 

Correlation - 




h 



wmax 



33.7 




d -°' 36 
P 



Remarks - None. 
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APPENDIX D: GENERAL EQUATIONS AND UNITS 



This appendix will be broken down into two sections. 

The first section will cover the general nondimensional heat 
transfer/fluids numbers used in this thesis. The second 
section will annotate the units to be used with each equation 
that is numbered within the main text. This requirement was 
necessitated by the use of different units (same variable 
nomenclature) among the various researchers cited in the 
thesis text. 

A. GENERAL NONDIMENSIONAL HEAT TRANSFER/FLUIDS NUMBERS 

1. Reynolds Number (Re) 

d u p- 
p o M f 

Re = 

y 

2 . Galileo Number (Ga) 

V Of ( p s - p f 1 9 

G 3. = ■ ■ ■■■ ■ 



Prandtl Number (Pr) 



c _ u 
pf M 



Pr = 
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4. Nusselt Number (Nu) 



h d 
w p 



Nu = 



Stanton Number (St) 



w 



St = 



c pf P f U o 



Biot Number (Bi) 



Bi = 



h r 
w sp 



sp 



B. EQUATION UNITS 

1. Equation (3.1) 



Ap e - lb/ft' 



- ft 

- dimensionless 

- lbm/f t 3 

- lbm/f t 3 



2 . Equation (3.2) 



'mf 



- dimensionless 

V - ft 3 
P 

Y f - lb/ft 3 

- ft 

- ft 2 

- lbm/f t 3 



J mf 
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p f - lbm/ft 3 
3 . Equation (3.3) 



G mf - lbm/ (hr-f t 2 ) 

D - in 

P 

y - centipoise 

P g - lbm/ft 3 

p f - lbm/ft 3 

Equation (3.4) 


D 


- in 


P 




D. 


- in 


1 




X i 


- dimensionless 


Equation (3.5) 


Re c 
mf 


- dimensionless 


G ^ 
mf 


- lbm/ (hr-f t 2 ) 


5 


- in 


P 




U 


- centipoise 


Equat 


ion (3.6) 


G r: 

mf 


- lbm/ (hr-f t 2 ) 


Equation (3.7) 


n 


- dimensionless 


G f 


- lbm/ (hr-f t 2 ) 


G 

e 


- lbm/ (hr-f t 2 ) 


Equation (3.8) 



r - dimensionless 
- dimensionless 



G f - lbm/ (hr-f t 2 ) 
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Ibm/ (hr-f t 2 ) 



9. 



10 . 



11 . 



12 . 



Equation (3.9) 



q • 


- 


BTU/hr 


h . T 
wiL 


- 


BTU/ (hr-ft 2 - °F) 


A 

m 


- 


in 2 


T 

w 


- 


°F 


T b 


- 


°F 


Equation (3.10) 


^wiL 


- 


BTU/ (hr-ft 2 - °F) 


k m 


- 


BTU/ (hr-ft- °F) 


p m 


- 


lbm/f t 3 


c 

m 


- 


BTU/ (lbm- °F) 


t 


- 


sec 


Equation (3.11) 


h w 


- 


BTU/ (hr-ft 2 - °F) 


k m 


- 


BTU/ (hr-ft- °F) 


°m 


- 


lbm/f t 3 


c m 


- 


BTU/ (lbm- °F) 


u s 


- 


ft/hr 


l h 


- 


ft 


Equation (3.12) 


h wL 


- 


BTU/ (hr-ft 2 - °F) 


k m 


- 


BTU/ (hr-ft- °F) 


c 

m 


- 


BTU/ (lbm- °F) 


s 


— 


1/sec 
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APPENDIX E: EQUIPMENT LISTING 



Device 


Manufacturer 


Model # 


Serial # 


Temperature 

Recorder 


Newport 


267B-TF2 


8500392 


Temperature 

Recorder 


Newport 


267B-TF2 


8500393 


Thermocouple 


Omega 

Engineering 


ANSI T 
Type S 
Sheathed 


— 


Water 
Manometer 
(30 in) 


Meriam 

Instrument 


Type W 
33KA35 


A69513 


Water 
Manometer 
(60 in) 


Meriam 

Instrument 


Type W 
A324 


44F554 


Flowmeter 


Fischer and 
Porter 


10A3565A 


6908A2165A1 


Variac 


General Radio 


Type V10 


— 


Air Blower 


Clements 

National 


G12 


114986 


Regulated 
Power Supply 


Lambda 

Electronics 


LK 345A FM 


C 1821 


Thermocouple 

Calibration 


Rosemount 

Engineering 


Bridge-9 20A 
Bath-9 13A 


110 


Strip 

Heaters 


Watlow 

Electric 


250 Watts 


8044 


Weight Scale 


Mettler 

Instrument 


Type H15 


119923 


Optical 

Microscope 


Bausch & Lomb 
Instrument 


31-32-45 


B-3539 RS 
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